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The description of the dynamics of correlated electrons in quantum impurity models is typically 
described within the nonequilibrium Green function formalism combined with a suitable approxi- 
mation. One common approach is based on the equation-of-motion technique often used to describe 
different regimes of the dynamic response. Here, we show that this approach may violate certain 
symmetry relations that must be fulfilled by the definition of the Green functions. These broken 
symmetries can lead to unphysical behavior. To circumvent this pathological shortcoming of the 
equation-of-motion approach we provide a scheme to restore basic symmetry relations. Illustrations 
are given for the Anderson and double Anderson impurity models. 



I. INTRODUCTION 

Describing the transport of electrons through an in- 
teracting region is a challenging task and typically in- 
volves the calculation of the dynamics of correlated elec- 
trons driven away from equilibrium 1 - - — . In general, this 
many body out-of-equilibrium problem cannot be solved 
exactly but for a few simple cases 4 - - — . Excluding recent 
developments based on brute-force approaches such as 
time-dependent numerical renormalization-group tech- 
niques^—, iterative 1 ^ - — or stochastic^ - — diagrammatic 
techniques to real time path integral formulations, wave 
function based approaches 1 ^, or reduced dynamic ap- 
proache o 20 ' 21 , all suitable to relatively simple model sys- 
tems, most theoretical treatments of quantum transport 
rely on approximations of some sort. One well studied 
approach is based on the nonequilibrium Green func- 
tion (NEGF) formalism otherwise known as the Keldysh 
NEGF or the Schwinger-Keldysh formalis m 22 ' 23 , which 
is widely used to describe transport phenomena^ 4 - - — . 

Based on the NEGF, an exact expression for the sta- 
tionary current through an interacting system coupled to 
large non-interacting metallic leads in terms of the sys- 
tem's Green function can be derived 2 ^: 

1 = £fj ds(Tr{f L (s-^ L )T L (e) 

x (G r (e)-G a (s))} + Tr{T L G<(e)}), (1) 
or equivalently 

I = ~ J deTr {£< (e) G> (e) £> (e) G< (e)} (2) 

where G r (G a ) is the retarded (advanced) Green function 
(GF) of the system, G< (G > ) is the lesser (greater) GF of 
the system, which will be defined later below. The lesser 
tunneling self-energy is given by = if L (e - fi L ) T L , 
where f k (e — fi k ) is the Fermi-Dirac distribution and T l 
is the matrix coupling the interacting system to the left 
reservoir with elements (T L ) mn = 2np k (e) t kn t* km (t krn 
is the hopping matrix elements between the system and 
the left reservoir). The calculation of the system's GF 
required to obtain the current (or other observables) is 



far from trivial, excluding simple noninteracting cases. 
Most applications are based on perturbative diagram- 
matic techniques to obtain G r , G a , G< and G > — . Al- 
ternatively, one can use the equation-of-motion (EOM) 
approach, which allows to deduce the system's GFs by 
deriving the corresponding equations of motion 2 ^ - — . In 
light of its simplicity, it has been used extensively to de- 
scribe transport phenomena such as the Coulomb block- 
ade^ 2 . and the Kondo effec t 30 ' 33 ' 34 , providing qualitative 
and in some cases quantitative results. When applied to 
interacting systems, the EOM for the GF gives rise to 
an infinite hierarchy of equations of higher-order GFs. A 
well-known approximation procedure is then to truncate 
this hierarchy, thus introducing a mean-field like descrip- 
tion to some observables. These equations for the GFs 
then need to be solved self-consistently for the resulting 
closed set of equations. Although successful, the EOM 
technique has its drawbacks^. 



In this paper we show that while a closure can always 
be obtained, it is not clear a priori whether it fulfills sym- 
metry relations that single particle GFs must obey. This 
failure can lead to solutions which are not physical, such 
as complex occupation of levels and even finite currents 
at zero bias. We also propose an approach to fix this defi- 
ciency by imposing a set of rules to reconstruct GFs that 
fulfill basic symmetry relations. Illustrations are given 
for the Anderson model 3 ^ at the Kondo regime and for 
the double Anderson model—. Our paper is organized as 
follows: in Sec. |TT] we describe the EOM approach and 
the single site and double site Anderson models. In Sec. 
IIIII we discuss symmetry relation for GFs and illustrate 
symmetry breaking for the aforementioned models with 
specific closures suitable to describe the Kondo effect. In 
Sec. HVI we provide a recipe to restore the basic symmetry 
relations within the EOM approach and discuss implica- 
tions for level occupancy and coherences, current, and 
sum rules for the Anderson model in the Kondo regime 
and the double Anderson model. Finally, in Sec. |V] we 
conclude. 
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II. EOM TECHNIQUE AND MODELS 
A. Equations of motion 

The EOM for the contour ordered GF— is obtained 
from the Heisenberg EOM for a Heisenberg operator 

f t A (t) = | [# (t) , A H (*)] + f i ff (t) , where in our 

case H (t) = H + V (t) . Here H stands for the one 
body noninteracting part of H (<), V (t) — H (t) — Hq, 
is the commutator. Let us consider a generic 



and 



A,B 



example. We define the contour ordered GF 



G (r 2 ,t 2 ,r 1 ,ti) = -- {Tc^h (r 2 ,t 2 ) (ri.ti) 



(3) 



where Tc is the contour time ordering operator and 
i s * ne system's annihilation (creation) field 
operator in the Heisenberg picture (in what follows we 
omit the H index). The EOM^ for G (r 2 , t 2 , ri, ti) can 
be written as (omitting the r dependence for brevity) 



92 (t a , «!)({*,**}) 
dtg 2 (*2,*) 



(4) 



«(t),V(t) 



where (ih-^ - e\ g 2 (t 2 ,ti) = 6 (h - t 2 ) , \A,b\ is 
the anti-commutator, and e is defined from the equa- 



tion, (t) = 



^ (t) , • For example, if H = 



and 'f = d, then e = e% — /i. Fol- 
lowing Langreth theorem 4 ^, we can change the contour 
integration in equation {2} to integration along the real 
time axis. This yields (see Sec. Mil for the definitions of 
the different real-time GFs) 



dtg r 2 (t 2 ,t)G r (t,t x ). 



G<(i 2 ,ii) = 5 2 < fe^i)({^,* t }) (5) 



dtg r 2 (t 2 ,t)G< (t.ti) 



dtg<{t 2 ,t) G a (Mi)- 



Co 



where G r (t 2 ,ti) is the retarded GF usually used to 
calculate the response of the system at time t 2 to 
an earlier perturbation of the system at time t\. 
G < (t 2 ,ti) is the lesser GF which plays the role of 
the single particle density matrix, and <&(t 2 ,ti) — 

* (t 2 ) , V (t 2 ) * t (ti)\ is a new GF generated 



by the EOM procedure. Depending on the Hamiltonian 



it can be a single particle GF or a many particle GF and 
can involve lead operators as well as system operators. 
In steady state, the GFs depend only on the difference in 
time, t = t 2 — t\, which is simpler to express in Fourier 
space 

G r (w) = gl (lj) ({*,**}) + .92 M G r M , (6) 



+g<(oj)G a (w). 



(7) 



To simplify the notation we denote the Fourier trans- 
form of G (t 2 — t\) — G (t) as G (w), i.e., functions with 
an argument "w" are Fourier transforms of their time- 
domain counterparts. At this stage one has to evaluate 
G (t 2 , ti) (G (t) in steady state). Except for very simple 
cases, where an exact closure can be obtained, writing 
the EOM for G(t 2 ,ti) will produce new and/or "higher 
order" GFs that need to be evaluated. This leads (in 
principle) to an infinite set of equations. The idea of the 
EOM method is therefore, to truncate this hierarchy of 
equations making a mean-field like approximation for the 
"higher-order" GFs through lower order functions. This 
is the Achilles heel of this method as there is no sys- 
tematic way to close the equations. Usually the approxi- 
mations have physical meaning within the regime of the 
problem at han d 33 i 41 ' 42 . In what follows we demonstrate 
that different approximations can sometimes break sym- 
metry relations that the GFs must fulfill. We will use two 
impurity models to demonstrate at what level of approx- 
imation the symmetry relations are violated and propose 
a scheme to restore symmetrization. 



B. The impurity models 

To illustrate the shortcomings of the EOM approach, 
we refer to the Anderson mode l 33 ' 36 ! 43 and the double 
Anderson model^l to represent two different degrees of 
complexity in correlated systems. As commonly used, 
we split the total Hamiltonian into three parts 2 ^: 



H = H 



sys 



Hba 



th 



Hi, 



(8) 



where H^ath describes the macroscopic leads (left and 
right contacts), H sys describes the system of interest 
(in our case the impurities), and Hi nt is the interac- 
tion Hamiltonian between the system and the leads. The 
contacts (leads) are modeled as infinite non-interacting 
fermionic baths^ 4 ^— with a Hamiltonian in second quan- 
tization given by 



Hbath — ^ £k,<rCk,a C k,<Ti 
cr.ke{L,R} 



(9) 



where e^.a is the energy of a free electron in the left (L) or 
right (R) lead, in momentum state k and spin a. The op- 
erator Ck,a (c\ is the annihilation (creation) operator 
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of such an electron. The form chosen for H sys depends on 
the system studied. For the Anderson impurity model 36 



o-e{t4} 



(10) 



Here n a = dJ a d a is the number operator of the spin a 
electron with energy e a and U is the repulsion energy 
between two electrons on the same site with opposite 
spins (intra-site repulsion). The second model we discuss 
is the double Anderson model^l 



H 



*tj* ■ ^ ^ £-111(7^111(7 ~b ^ u rn Tt rn -^7i n i^ (n) 

o-,me{a.,fi} m 



where the first two terms on the R.H.S are similar to 
the Anderson impurity model Hamiltonian (extended to 
2 sites), V*2 is the repulsion energy between two elec- 
trons on different sites (inter-site repulsion), and h" a ^ is 
the coupling strength for electron hopping between the 
two sites. The interaction between the system and the 
contacts is simply given by the tunneling Hamiltonian 47 



Hi, 



(12) 



The parameter t% m represents the coupling strength be- 
tween the system and the leads, and the index m runs 
over the site index {a, /3} in the double Anderson model. 



The retarded GF can be used to calculate the response 
of the system at time t 2 to an earlier perturbation of the 
system at time t\ and is proportional to the local den- 
sity of states, while the lesser GF is also known as the 
particle propagator and plays the role of the single parti- 
cle density matrix. From equation ([I]) it is obvious that 
in order to calculate the stationary current the retarded, 
advanced and lesser GFs are needed, thus, the current is 
expressed in terms of the local density of states and the 
occupation of the system. Using the given definitions it 
is clear that the following relations must hold: 

Gtffati) = (ii,t 2 ))*, (14) 

G a p (*2,*l) - G%p (*2,tl) = G aj3 (*2,fl) — G a p (t 2 ,t%) . 

In steady state these relations can be rewritten in Fourier 
space as: 



G r a f) (<*>) = (Gp a (u)))* , 
G^H = -(G^H)*, (15) 
Glp M - G a p M = G> ap (w) - G% (lj) . 



In what follows we show that these relations do not hold 
when the GFs are obtained by the EOM technique with 
an arbitrary closure. 



B. The Anderson model 



III. SYMMETRY BREAKING IN IMPURITY 
MODELS 

A. Definitions and symmetry relations 

In the Keldysh formalism the two time NEGF is de- 
fined on a contour. In accordance with where on the 
contour the two times are placed one can define six real- 
time GFs^; the time-ordered G', anti-time ordered G 4 , 
lesser G < , greater G > , retarded G r , and advanced G a : 



Glpfah) = -^(i 2 -*i)^* a (ia)*J(ti) 



+~6(t 1 -t 2 ){*} (h)y a (t 2 )), 



Gap(h,h) = -^(ii-* 2 )(* a (ia)*J(ti) 



I 
i 

t 



+^(ia-*i)(*J (*i)* Q (t 2 )>, 



G aP (*2,*l) 



-^t { t l) i, a ( t2 )), 



(13) 



G a p{t 2 M) = -^(*2-*l)({^a(*2),^(*l)}}, 

G a p(t 2 ,h) = ^(ti-t 2 )({* a (t 2 ),^(fi)}' 



Following the derivation in Refs. 28 33,47 we define the 
following contour ordered GF: 

G aa (t,t') = - l -(T c d a {t)dl{t')), (16) 

G 2 (t, t') = ~ (T C n s (f) d{j (i) 4 (;')) , (17) 

where a is the opposite spin of a. Various approximate 
decoupling procedures can be applied to the many par- 
ticle GF— . Here we follow the approximation scheme 
used in Refs. [2811331 where all electronic correlations con- 
taining at most one lead operator, are not decoupled 
and their EOM are calculated. Higher order GFs in- 
volving (opposite) spin correlations in the leads are set 
to zero, and the remaining higher order GFs involving 
lead and system degrees of freedom are decoupled such 

that F 2 (£, t') = -i (T c c{ nS (t) d a (t) c qms (t) 4 (f )) = 

-hqSmnfk (e n ~ Mfe) G<j<j (*,*')■ Tne resulting EOMs (in 
Fourier space) are: 

(fiw - e a - E (w)) Gaa (w) = 1 + UG 2 (w) , (18) 

G 2 H = (^-e ff -[/-S H-S 3 ( W r 1 

x «n ff ) - Ei (w) G CTCT (w)) , (19) 
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where (n s ) = 



i ke{L r} hu)~e k — ls exac t self-energy for the noii- 



interacting case, £1 (w) and S3 (oj) are the self-energies 
due to the tunneling of the a electron , and are given by 



fce{i,-R} 



1 



fej + e k . s - e a - e s - U 
1 



(20) 



f^J - £k,a - Ecr + S d 



J =1,3 



with = / fc (£ fei<7 - Hk), A k °> = 1, and f k (e k .a - Mfc) 
is the Fermi Dirac distribution. To show that these 
set of equations break the symmetry relation G^ a (cj) = 
— (G^ct (w))* we define 



,(3) 



S 4 (cj) = E M + S3 (w) 



E (w) ' 



52 (w) 



— e a — U — £ 4 (cj) 



(21) 



(22) 



(23) 



With these definitions equations (Tig)) and (fT9")) can be 
rewritten (omitting (u>) for brevity) as: 



G 



gUG 2 , 



G2 = 92 {ns) - g2^iG acy . 



(24) 



(25) 



Substituting equation (|25|) in equation (|24[) and applying 
the Langreth rules we find that the lesser GF is given by 

G< ff = <?< + .9 r f/P r 5 2 < (n»> + 3 < ^ Q 52 (K) - £?ff Q ) 
-9 r UP r g r 2 (E r l9 < + E< ff Q ) - g r UP r g<Z a ig a 

-g r UP r g^\g<UP a g a 2 ((n ff ) + E? 3 Q ) (26) 
-g'UPrfiVfg'UP'gS (K) + £? 5 a ) 
-g r UP r g?X a ig a UP a g$ {(n s ) + £? 5 a ) , 



where P' r,a = 



.9 



g r T,<g a , and 3^ 



32^4 52 ■ Applying the principle of reductio ad absurdum 
we assume G^ a is imaginary. Since it must hold for any 
real value of (n s ) between and 1, we argue that the 
term 

A 1 = g r UP r g<(n s )+g<UP a g%(n s ) 
-g r UP r g r 2 Y.{g<UP a g a 2 {n s ) 
-g r UP r g r 2 iZ<g a UP a g a 2 {n s ) (27) 
-g r UP r g<^g a UP a g^ (n s ) , 

must be imaginary. Moreover, Since A\ must be imagi- 
nary for any value of U the term 



should be imaginary as well. Using the fact that U and 
(n s ) are real quantities and by definition g 2 and g < are 
imaginary, for A 2 to be imaginary the following must 
hold: 



(29) 



or in other words, we demand that 5ft (A 2 ) — 0. One can 
then show (see online supporting material for more infor- 
mation) that, in fact, the equality in equation (|29p does 
not hold, namely, G< a (w) is not an imaginary function 
and the relation G< a (us) = — (G< CT (u>))* is not satisfied. 
In turn, this implies that (n a ) (the occupation number) 
is a complex number, which of course is not physical. Fol- 
lowing the same derivation one can show that G„ a (u) is 
not an imaginary function either. All the other relations 
given in equation (|15p are fulfilled. 

If one is only interested in the Coulomb blockade 
regime, it is not necessary to go to the level of approx- 
imation presented here (which is essential to obtain the 
Kondo effect). For the Coulomb blockade regime one can 
turn to the approximation presented in Refs. 32, where 
on top of the approximations described above we also ne- 
glect the simultaneous hopping of electron pairs to and 
from the system. This approximation does not violet 
the symmetry relations of the single particle GF (see on- 
line supporting information for further discussion), but 
as pointed above, it does not reproduce the Kondo peaks 
at low temperatures. 



C. The double Anderson model 

For the double Anderson model we follow the deriva- 
tion given in Ref. [5(J and define the following contour 
ordered GF 



G 



Vp (*> = ~\ (Tcd aa (t) d\ a (t')) , (30) 



1 Tcrcr 



(t, t') = -- (T c n aT (t) dp a (t) dt (f)) , (3i) 



where r = a, a . The approximations used in 
Ref. HO are: (a) neglect the simultaneous hopping 
of electron pairs to and from the system, (b) as- 
sume that F 2 (t,f) = -l(T c c H<J {t)n{t)d\ a {if)\ « 

-I E 1=a , t%a I dtig k (t, h) (r c d ia (to n (to d\ a {t')) 

where n (t) is the number operator of one of the 
electrons of the system, and (i/ijj — £fc CT ) 5% (t, t0 — 
5[t — ti), and (c) higher order GFs of the form 



h 

pled 



are decou- 



to 



n ll<7 (t) ng, r (t) d a><7 (t) d'p t(r (0) 

(t))(r c n s , T (t)d a , a (t)dl J0)) - 



An 



g r UP r g^ {n s ) + g < UP a g 2 l (ng) , (28) i ( n s,* (*)> (^,^ (*) <W (*) dj.tr (°)) ■ These a PP roxi - 
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mations lead to the following results 

G% (w) = (hw - e a , a ~ So Mr 1 x (6% 

+hZpGffi (w) + t/ a G3 (w) (32) 

+cjg^ («) + c/g^ ( W )) , 



G^ (w) 
G^ (w) 
G^ H 

G^? (w) 



W/3 



(fej - e Qff - C/„ - V a 7 (n^) - V£| («/?») ~ So (a;)) -1 

x [hlpGffi («) + (n aff ) (T^fG^ H + ^TG|S? H)] , 

(fiw - £(3ct - (n fjs ) - Vft (rw) - Vfc - So M) -1 

x [(n a9 ) + h$ a G%? p (w) + («a ff ) (t^Gf^ (w) + P^G^ (w))] , 

(fiw - Bj9<F - Up ( nfl9 ) - Vf° (n aS ) - - S (w))" 1 

x [<tw> + M + <»W> (Wj^ (w) + H)] , 

(Tiw - £ Q(T - f7 a (n aff ) - V£p (np a ) - V°p - S (w)) 1 

x 'Kffifi (w) + (r^) (t/aG^ M + V^G^ M)] , 

(fiw - e Qff - (7 Q (n as ) - V£j - V% ~ So (a;))" 1 

x - + ^G^ (w) + (n^) (V Q G^ (w) + V£JG^ 

(fe - £(3ff - Up - Vf° (n aff ) - VJ a CT (rw) - S (u;)) -1 

<*/>*> + ^ a G|S? H + M (VpfGZpf) (w) + t^fC^ (a,)) 



(w) 



(33) 



We now show that given this set of equations, the sym- 
metry relation (g^ (wH =(fe(u)J ] is not sat- 
isfied. By applying the Langreth rules we can find the 
retarded and advanced projections of the single particle 
GF (equation Q32p). For simplicity we derived them for 

the case where VfJ — 0. Define 

1 j 



Given these definitions, the retarded and advanced GFs 
are given by: 



{9iY 



hui — Ei.n — S n 



huj - Ej, a - Uj (nj,a) ~ S '' 



(34) 
(35) 
(36) 



{G7pY = (l-(g a ) r h^(gp) r hl a ) 1 (37) 
x (( 9a ) r hip (gp) r + (g a ) r Kp (gp) r Up {&$%) r + (g a ) r U a (G^) r ) , 
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x (1 - {g%) r h° ia {gf a ) r hip)' 1 {g%) r h°p a {g$l) r (n ps ) U a ) ~ l (38) 

x ((C) r Kp (gZpY (««f) + (gZY Kp (9apY to*) U P 
x (1 - {g$l) r h°p a {g$l) r hip)' 1 {g%) r (np- a )) , 

{G$^) r = (1 - {g%°p) r h°p a {g$l) r Kp - {g%%) r h% a (g$l) r (n [jg ) U a 

x (l-te) r ^(^) r ^a) _1 te) r ^(^) r (^)^) _1 (39) 
x {{MY to*) + iMY h p a {MY to*) u a 

x (i - (gZY Kp {9%Y h % a Y x (gZY Kp {MY to*)) , 

(g%y = (i-(gp) a h° Pa ( 3a ) a hz fl y 1 (40) 

x {(gp) a h% a (g a ) a + (gp) a h" Pa (g a ) a U a (G™)" + (gpf Up (G$ZY) > 

{ G PPaY = (l - (gppY h Pa (9/3aYKf3 - (gppYhpa (gpaY to*) u a 

x (l-^) ^^) ^)" 1 ^) ^^)"^)^)" 1 (41) 

X ((.9f?r ft? Q (C) a + (flf?)" (4a) " <»*r> ^ 

x (i - {gZYKp (sS?)°^a) _1 (s£) a - 

(^Z) a = (l-(^) a ^(5^ a ^a-(flaa)°^(5SS) a <»a ff )^ 

x (l-(^) a ^ Q (4S) a ^) _1 (4^) a ^ Q (^) a (^)^) _1 (42) 

x ((<&)" W + 0C)° hip (gZpY to*) Up 

x (1 - {g%) a h°p a {gf a ) a hip)' 1 {g%) a h° pa (g$lY (np*)) . 



Substituting the equations for y^ZxP ( W )J ' 

(G|^(w)) r , (G|^(w)) a and (G*£ (u)) a into 

equations (|37|) and (|40p. respectively, and com- 
paring the resulting expressions we find that 

(G^H) r ^ ((G^H) Q )* (see online support- 

ing material for more details). Moreover, we find that 
none of the symmetry relations in equation (|15l) hold. 
In the following section we propose a symmetrization 
scheme that restores all the symmetries of the single 
particle GF. 



IV. SYMMETRY RESTORATION 

A. Guidelines to restore symmetry 

The customary route to calculate the NEGF is as 
follows: (a) calculate the retarded GF and use it to 
obtain the advanced GF (by demanding G% (ui) = 

{^G r p a (w)^ ). (b) Calculate the lesser/greater GF 

and symmetrize the lesser/greater to fulfill the quan- 
tum Onsager relations^, hence obeying G^'p" (ui) — 

— (Gp'a (u>)) . In most applications of NEGF the ad- 
vanced GF is not directly calculated and thus, the sym- 
metry breakage does not always stand out. In fact, this 
common procedure restores the relation between the ad- 
vanced and retarded GF and between the lesser/greater 
and their complex conjugate, but does not necessarily re- 
store the relation G r a p (u>)-G a a p (ui) = G>p (u))-G<p (ui). 
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It can be shown that violation of the latter leads to 
violation of the fluctuation dissipation relation, G < = 
—feq (s — l^eq) (G r — G a ), at equilibrium. This oversim- 
plified procedure can result in different values for the 
currents depending on how it is calculated, cf. equation 
|[TJ or equation j2]). It may also lead to finite currents 
at zero-bias voltage (see Sec. II V CI for more), which is 
physically incorrect. 

In order to restore the symmetry relations that are 
imposed by the definitions of the GF (cf. equation ITS]) , 
we suggest the following procedure: 

1. Calculate the retarded/advanced GFs matrices 
(G'/G a ) separately and use them to define 
"new" retarded/advanced GFs matrices G r = 
\ (G r + (G a ) f ) and G a = ± (G a + (G r ) f ) = 

2. Use the "new" retarded/advanced GFs matrices 
^G r /G a ^ to calculate the lesser/greater GFs 
matrices (G < /G > ). Again, use them to de- 
fine "new" lesser /greater GFs matrices G <:> = 
l(G<<>-(G<.>)t). 

3. Calculate the two anti-Hermitian matrices A = 
G> - G< and B = G r — G a . Define the difference 
anti-Hermitian matrix C = A — B, and redefine 
the retarded and advanced GFs G r = G r + % , and 
G a = G a - f . 

The resulting GFs (G r , G a , G< and G>) obey all sym- 
metry relations of equation (p~5|) by construction. Note 
that if the original GFs obeyed the symmetry relations to 
begin with, our symmetrization procedure will not alter 
them in any way. 

We now turn to perform detailed calculations for both 
the Anderson and double Anderson models. For the An- 
derson model, we use the closure described in Sec. IIII Bl 
while for the double Anderson model we use the closure 
described in Sec. IIII CI The resulting EOMs were solved 
self-consistently in Fourier space with a frequency dis- 
cretization of — 2 14 — 2 16 depending on the model pa- 
rameters. Typically, < 15 self-consistent iterations were 
needed to converge the results. Convergence was declared 
when the population values at subsequent iteration steps 
did not change within a predefined tolerance value cho- 
sen as 10~ 6 . For each set of calculations we have applied 
the above symmetrization scheme and compared the re- 
sults to those obtained without restoring symmetry, as 
detailed for each model. 



B. Anderson impurity model 

First, we address the effects of symmetry breakage in 
the Anderson model. The closure used is sufficient to de- 
scribe the appearance of the Kondo resonances at low 



temperatures, as seen in the upper panel of figure [TJ 
where we plot the density of states as a function of en- 
ergy for several temperatures, all calculated with symme- 
try restoration. The development of Kondo peaks in the 
density of states as the temperature decreases is clearly 
evident, signifying a regime of strong correlations which 
is qualitatively captured by the simple EOM approach 
when symmetry is restored. 

In the lower panel of figure [TJ we show one of the main 
flaws of the EOM approach for the Anderson impurity 
model, where we plot the value of (n-j-) as a function of 
the source drain bias voltage with and without symmetry 
restoration. The most notable effect is the appearance of 
an imaginary portion to (n-j-) as the source drain bias 
voltage is increased. To obtain the results, without sym- 
metry restoration, only the real part of (n CT ) was used to 
converge the self-consistent equations for the GFs. By ap- 
plying the symmetrization scheme proposed in Sec. IIV Al 
to the lesser GF calculated in Sec. IIII Bt we restore the 

relation G^'^ (cj) = — (g^^ {<-^)j ■ This is sufficient to 

obtain a real value for (n^), as clearly shown in the lower 
panel of figure [TJ All other symmetry relation are not vi- 
olated here and thus, our symmetrization procedure does 
not affect them at all. Interestingly, taking only the real 
part of {n a ) provides identical results when compared to 
the results obtained after the full symmetrization pro- 
cedure. However, this is only true for the simple case of 
the single site impurity model and does not hold for more 
complex systems. 



C. The double Anderson model 

We now turn to discuss the impact of symmetry break- 
ing for the double Anderson model. This system is more 
involved compared to the single site Anderson model and 
thus, the level of closure used is somewhat simpler, as 
explained in Sec. IIII CI While for the case of a sin- 
gle site Anderson model only the relation G„i (oj) = 

— [G*p'a (<^)) breaks down, in the double Anderson 
model we find that all 3 symmetries described by equa- 
tion [15] are violated. This can be traced to the more 
complex form of the Hamiltonian for the double Ander- 
son model, where each site is only coupled to one of the 
leads and transport in enabled by the direct hopping term 
between the two sites. 

Similar to the case of the Anderson model, as a re- 
sult of symmetry breaking the occupation of the levels 
(n Q(T ) is a complex number. In addition, the coherences, 

PZ = -£f-oo ( G *p(")) <d "> should also fulfill cer- 
tain symmetry relations, such as p^ = ■ In fig- 

ure [5] we plot the real and imaginary parts of p^Z and 
Pfia f° r t ne case where the symmetry procedure has been 
applied (left panels) and for the bare case (right pan- 
els). The upper panels show the imaginary part of p^l 
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Figure 1: Upper panel: Density of states in the Kondo regime 
for nonequilibrium situation of the spin up electron after sym- 
metrization for different temperatures. Parameters used are 
similar to those used in Refs. 34,52 (in units of T = Fl + Tr): 
fj,L = 3/10, /ijj = 0, Ei,f — — 2, and U = 10. The bands are 
modeled as a Lorentzian with a half bandwidth 100. Lower 
panel: Occupation of the spin up electron before ("bare") and 
after ("sym") symmetrization. As can be clearly seen the 
real part of the observable (n-f) is not affected by the sym- 
metrization, and only the non physical imaginary part disap- 
pears. Parameters used (in units of U): Tl^ = Tr^ = 0.3, 
T i4 = T R:l = 0.05, e t = 0.2, £ir = -0.2, /3 = 4 and U = 1. 



and p a g a , which should show a mirror reflection about the 
zero axis (shown as thin solid line). This is, indeed, the 
case when symmetry is restored, however, it is destroyed 
when symmetry breaks down, in particular as the source 
drain bias increases. A more dramatic effect is shown 
for the real part of p°p and p a ^ a (lower panels). The 

two curves representing 5ft {ji'^j and 5ft {p'p'C] snou ld be 
identical (left panel when symmetry is restored) but are 
quite distinct when symmetry is not obeyed (right panel) . 

In figure [3] we plot the current as a function of the 
source drain bias voltage for the double Anderson model. 
The current can be obtained from equation ([!} (dashed 



Figure 2: The imaginary (upper panels) and real (lower pan- 
els) parts of = (d^dairj (dashed line) and ppl = 

(f$ aa d{icr) (solid line) calculated before (right panels) and af- 
ter (left panels) symmetry was restored. The solid thin line 
in the upper panels marks the zero axis. As expected, after 
symmetry restoration (left panels), 3 (Pap) = {Ppa) an ^ 
%t (.Pip) = ^ (p < Ba)i while before symmetry restoration (right 
panels) these equalities are violated. Parameters used for the 



simulations in units of U = U a — Up are: — ± a ^ 
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0.25, Vgl = 0.1, e at = e al = 0.1, e 0t = e n = -0.175 and 
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line) or from equation @ (dotted curve). In the limit 
of infinite hierarchy in the EOM approach the two for- 
mulas should coincide. However, when approximations 
are introduced or when the hierarchy is truncated, the 
calculation of the current based on the two different 
formulas will coincide only if the symmetry relation 

M) r - {Gil H) Q = (g% (^)) > - (g% {^Y 

is preserved. Indeed, in the case of a single site Ander- 
son model, even if symmetry is not restored, this relation 
holds and the two calculations yield identical values for 
the current. However, in the present case, all 3 symme- 
try relations are broken and thus, equations ((I) and @ 
give different results for the current, as clearly evident in 
figure |31 More significantly is the fact that equation ([T]) 
produces a finite value for the current even when the bias 
is zero, indicating the break down of the fluctuation dissi- 
pation relation. When symmetry is restored (solid curve) 
the two calculations are identical, as they should be, and 
the violation of the fluctuation dissipation relation is also 
resolved. 

The symmetrization scheme proposed here is not a 
"magic cure" and, in fact, does not resolve all issues of 
mater. It is well known that the lesser and greater GFs 
should obey a simple sum rule where the integral over 
the difference of their diagonal elements should always 
sum to 1: 



S «° = ^f d£ (( G - W) < ~ ( G ™« (*))*) = 1- 



(43) 




Figure 3: I-V curves calculated using equations {J) and © 
before (dashed and dotted lines) and after (solid line) ap- 
plying the symmetry procedure suggested in Sec. I IV Al As 
can be clearly seen, before symmetrization, calculating the 
current via the two different but equivalent formulas provide 
different results, one of which is not physical (dashed line) 
as the current is finite for Vsd = 0. The latter result sug- 
gests that the "unsymmetrized" GFs obtained through the 
EOM disobey the fluctuation dissipation relation. Parame- 
ters used for the simulations in units of U — U a = Up are: 

-pL t^R taR r) f\(\r)r pi? pi? pL 

1 at — 1 c4 — 1 0t — 1 Pi ~ U.UU/O, 1 a -f — 1 a ^ — 1 p-f — 

Tk = °> h *s = Kp = 0.25, VZ T p = 0.1, Ea t = e al = 0.1, 
ept = £pi = -0.175 and /3 = 80. 



In figure U we plot the sum rule as given by equation 
([4"3")l for the double Anderson model where symmetry has 
been restored. A similar plot for the single site Anderson 
model yields a value of 1 regardless of whether symme- 
try has been restored or not within the closure discussed 
above. However, in the case of the more evolved dou- 
ble Anderson model, even when symmetry is restored 
and the GFs obey all 3 relations described in equation 
(|T5]l . the sum rule is violated. Nonetheless, the sum 
Yl^cM? — N e , where N e is the total number of electrons 

in the system at maximal occupancy, is indeed preserved 
when symmetrization is restored. 



V. SUMMARY 

In this paper we have addressed the problem of sym- 
metry breaking and restoring in the EOM approach to 
NEGF formalism. This formalism is based on deriving a 
hierarchy of equations of motion for the system's Green 
functions and truncating this hierarchy at a desired (or 
tractable) order. Despite the uncontrolled approxima- 
tion introduced by an arbitrary truncation, the closed 
set of equations is often used to describe the complex 
dynamics of correlated systems, including the Coulomb 
blockade and Kondo effect. 

One shortcoming of the EOM approach, which has 



1.1 



0.1 0.2 0.3 0.4 

V SD /U 

Figure 4: £w (dashed line) and Sp a (dotted line) calculated 
from the "symmetrized" lesser and greater GFs as a function 
of the source drain bias voltage. The exact result should have 
been 1 (as marked by the solid line) . Parameters used for the 
simulations in units of U = U a — Up are: = = Ff^ = 

r« = 0.0025, r« = r« = r£ t = = o, k p = h' aP = 

0.25, V$ = 0.1, e at = e al = 0.1, ep t = e pi = -0.175 and 
13 = 80. 



been the focus of the present study, is the fact that, a 
priori, for most situations it is impossible to determine 
whether the solution of the closed set of equations satis- 
fies symmetry relation between the retarded, advanced, 
lesser and greater Green functions imposed by definition. 
For example, we have shown that for the Anderson model 

the relation (u) = — [Gp^ (w)J breaks down for 

a closure that is often used to describe the dynamics 
near the Kondo regime. We have also demonstrated that 
for the double Anderson model all 3 symmetry relations 
given by equation (|15p break down for a lower level of 
closure. This faulty of the EOM approach leads to un- 
physical behavior such as complex level occupations and 
finite current at zero source drain bias (depending on how 
the current is evaluated). 

We have also proposed a procedure to circumvent 
this deficiency by imposing symmetrization to the Green 
functions in such a way that all 3 symmetry relations are 
restored. The strength of the proposed approach is that it 
does not alter the GFs if symmetry is not broken. While 
this procedure eliminates some problems of physical im- 
portance and leads to real level occupations and vanish- 
ing current at zero source drain bias (irrespective of how 
the current is evaluated), certain sum rules are still vio- 
lated, indicating other problems with the EOM approach. 
Nonetheless, the symmetrized version of the EOM tech- 
nique still describes the appearance of the Kondo peak 
and, as will be shown in future publication provides a 
quantitative description of the resonant transport for the 
double Anderson model even in the strong inter-dot cou- 
pling limit. 
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Supporting Information 

VII. FULL DERIVATION OF THE BROKEN SYMMETRY IN THE ANDERSON MODEL 

Here we present in greater detail the breakage of the relation G< CT (uj) = — (G^ a (uj)) described in subsection The 
Anderson model in the manuscript. We demonstrate that G< CT (uj) is not an imaginary function. We start by 
defining the following contour ordered GF: 

G aa (t,t') = -^(TcdAt)di(t')), (44) 

G 2 (t,f) = -^(T c n a (t)d a (t)d a (t')), (45) 

where a is the opposite spin of a. The resulting EOMs (in Fourier space) under the approximation scheme discuss in 
the manuscript are: 

(foj - e a - S (w)) G aa (w) = 1 + UG 2 (uj) , (46) 

G 2 (uj) = (^-£ ff -C/-EoH-S 3 H)~ 1 ((7i ff )-SiHG CT(T H), (47) 
We define the following GFs and self-energies: 

9^) = ~ h 1 v t V ( 48 ) 

Tiui - e a - S (uj) 

92 («) = 1 , (49) 

hid — e„ — U — S4 (uj) 

S H= h '^' 2 ■ ( 5 °) 



i,fee{L,_R} 



(51) 

S 4 (w) - So («) + S 3 (w) . (52) 

with = /fc (ej fc,<r ~ Mfc)> -^jt = 1) an d /fc ( £ i,fc,<r ~ Mfe) i s the Fermi Dirac distribution. 
Rewriting the equations of the GFs in terms of the above definitions gives: 

G 2 (uj) = 52 M (n-sr) - .92 (w) Si (uj) G a( j (w) , (53) 

Gcrcr (uj) = g (uj) + g (uj) UG 2 (uj) , (54) 
We then merge equations (|54|) and (|53|) to get: 

G 2 (uj) = g 2 (uj) (n s ) - g 2 (uj) Si (uj) g (uj) - g 2 (uj) S x (uj) g (uj) UG 2 (uj) . (55) 
The advanced GF can be extracted from the contour ordered one (equation (155)0 by setting^ u) — > uj — ?'0 + 



G a 2 (uj) = g a 2 (uj) (n s ) - g a 2 (uj) E« (uj) g a (uj) - g a 2 (uj) E« («) 5 Q (w) t/G° («) , 



(56) 
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For brevity, we omit (w) and rewrite equation (|56p as: 

G a 2 = (l+g^ a 1 g a U)- 1 g^(n^-(l + g^g a U)- 1 glE a 1 9 a . (57) 

Define: 

P ra = (1 + g^^g^uy 1 , (58) 

G a 2 = P a g c 2 l (n 9 ) - P a g^ a l9 a , (59) 
Using Langreth theorem, the lesser projection of G 2 can be evaluated: 

G 2 < - ,9 2 < (n 9 ) - ( ff2 Si.9) < - ( ff2 E l5 ) r C/G< - ( 52 E l3 )< UG a 2 , (60) 



with 



,9< = gH<g*, (61) 



,9< = <££^, (62) 

and the lesser self energies are defined as in Ref. 

x< = x< L + x< R = i(f L r xL + f R r xR ), (63) 



where 



r x , = -2/m(E^), (64) 



and E^f stands for the retarded ("r") or advanced ("a") self-energies. Substituting equations (f6"Tj) and (f6"2")) into 
equation ([60)1 , the lesser projection of equation ([55]) is given by: 

G< = P r g< (n 9 ) - P r g r 2 V[g< - P r g r 2 E<g a - P^E?^ 

-P r g r 2 Z r ig <UP a g2 > (n 9 ) - P r g r 2 Z<g a UP a g2 > (n ff ) - P r 9 <E a l9 a UP a g^ (n 5 > 
+P r <££J S < E/P a <£E? S ° + P r g r 2 i:<g a UP a g^g a + P r g<^ a ig a UP a 9 ^ a ig a , (65) 

The lesser projection of equation (|54[) can now be written as 

G< a = g< + g r UG< + g<UGl (66) 

Using our results for G 2 and G^ we find 

G< CT = ,9< + g r UP r g< (n 9 ) + ^UP^ ((n s ) - E? ff a ) - g r UP r g<^g a 

-g r UP r g r 2 (Z r l9 < + E< 5 Q ) - g r UP r g r 2 ^g<UP a g a 2 {(n 9 ) + £? ff a ) (67) 
-g r UP r g r 2 i:<g a UP a gZ ((n 9 ) + £? 5 Q ) - g r UP> ' g<^g a U P a g a 2 {(n s ) + E?.g a ) 

Applying the principle of reductio ad absurdum, we assume G^ CT is imaginary. Since it must hold for any real value 
of (n s ) between and 1, we argue that the term 

A l = g r UP r g< (n s ) + g<UP a g a 2 (n s ) - g r U P r g r 2 V l9 <U P a g a 2 (n s ) 

-g r UP r g r 2 E<g a UP a g^ (n 9 ) - g r UP r g<E a l9 a UP a g^ (n 9 ) , (68) 

is imaginary by itself. Moreover, Since A\ must be imaginary for any value of U, the term 

A 2 = g r UP r g< (n 9 ) + g<UP a g% (n 9 ) , (69) 

should be imaginary as well. Using the fact that U and (rig) are real quantities and by definition g 2 and g < are 
imaginary, for A 2 to be imaginary, one requires that its real part vanishes, i.e.,: 

Im (g r P r ) g< + Im {P a g a 2 ) g< = 0. (70) 
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In other words the equality 

Im (g r P r ) g%Z<fi = -Im (P a g a 2 ) g r ^g a , (71) 

must hold for the assumption that G^ a is imaginary to be satisfied. Using the definitions for g and 172 the last equality 
can be rewritten as: 

Im ( Q r P n -if(u,)Im(£ r 4 ) _ Im(pa a) »/»7m(E£) 

{hw-e^-Uf + (7m(£^)) 2 y2 ' (^-£ ) 2 + (7m(£ ')) 2 ' 

where eo = e a + Re (Eg) and £4 = e a + Re (£4). Starting with the L.H.S. of equation (|72p. we look at g r P r 

9 r 1 



g rpr 



1 + gVZWU huj-e - E r + 

?kj — e a — U — £4 
(huj -e a - Eg) (huj-e a -U- SJ) + £^J7 

fru — E4 — U — i (JmSJ) 



(/kj - e - i (/mEg)) (TtD - s A - U - i (JmEJ)) + U Re (££) + 17 ■ i (7m££) ' 
Denote ao = fiw — £0, ai = Tfe (££), 04 = fta; — £4 — U and 6a; = Jm (££) with x = 0, 1, 4 

04 — 264 



(ao — i&o) ( a 4 — ^4) + C^ a i + iUb\ 
04 — 264 



aoa4 — 6064 + C^ai — i (^0^4 + 0460 — Ub\) 
Denote D = a$a^ — b^b^ + Ua\ and E — ao&4 + 0460 — Ub\, so we can rewrite equation (|74l) as 



Finally 



The L.H.S. of equation ([72]) is thus 



^.92 



1 + g^1g a U nu - £ CT - U - E? 



S?(7 



fiw — £ CT — Eg 



(73) 



(74) 



r r 04 — 164 (04 — 164) (7? + iE) a^D + b±E — i (b±D — a^E) 
9 = D - iE = D 2 + E 2 = WT~E 2 ' 1 ' 



*™(9 r n = - b -§^ (76) 



Im(g r n Z ifMpm = tjph 2 hD - a f . (77) 

(fiw - £4 - C/) 2 + (Im (£4)) (a 4 ) + (64) ^ + ^ 
Now we turn to analyze the R.H.S. of equation (f72"j) . We start with evaluating P .^ 



(78) 



{huj-e a -U- £2) (fiw - e a - Eg) + Ef C/ 

ao + ^0 

(ao + «6o) (d4 + 164) + J7ai — zC/6i 
_ a + ib 

0004 — 6064 + Ua% + i (ao&4 + (1460 — Ub%) 
_ ao + i&o _ («o + «M (-D — iE) _ a D + b E + i (b D - a Q E) 

D + iE ~ D 2 + E 2 ~~ D 2 + E 2 ' 

To go from the second line to the third line in equation (|78[) we used Im (££) = —Im (£"). Finally: 

Im ™)=*§^£. (79) 
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The R.H.S. of equation (JJ3) is thus, 



(r^- £o ) 2 + (/m(E 5 )) 2 £ 2 + £ 2 (a ) 2 + (6 ) 2 



The equality fequation (|T2"|0 now reads: 

64 b 4 D — a 4 E b$D — atjE bo 



(a 4 ) 2 + (6 4 ) 2 ^ 2 + ^ 2 ^ 2 + ^ 2 M 2 + (M 2 ' 

or 

b\D - a 4 b 4 E _ b%D - a Q boE 
(a 4 ) 2 + (6 4 ) 2 " (a ) 2 + (bof ' 



(81) 



(82) 



Substituting D = aoa 4 — bob 4 + t/ai and E 1 = ao& 4 + a 4 &o — Ub\ one can easily show that the equality does not hold. 
Thus, G< CT (w) is not an imaginary function as it should be by definition. 

In the paper we argued that a simpler closure (as used for example, in Ref. I32I ) will not violate the symmetries of 
the GFs. In what follows we show that under the simpler closure, indeed, G<„ (uj) = — (G< a (to))* . Our starting point 
is the same. Define the contour ordered GFs: 

G (T(T (t,t') = ~{T c d (T (t)4(t')), (83) 

G 2 (t, = -~ {T c n s {t) d(j (t) 4 (f)) , (84) 

Following the approximations of Ref. |32|, the resulting EOMs (in Fourier space) are: 

{tyjj-e a - S H) G ff(T (w) = 1 + UG 2 (w) , (85) 

G 2 (w) = (fc-^-CZ-EoHr 1 ^). (86) 

We define 

■9H = i 1 v ( v ( 8? ) 



52 H 



— £cr — U — So (w) ' 

Rewriting equations ([55)1 and (|8"6")l in terms of the given definitions we get: 

G 2 (w) = 92 M (n*) , (90) 

G aa (w) = .9 («) + g (w) C/G 2 (w) , (91) 
Substitute equation (|9"0")) into equation (|9"Tj) 

G M = 3 (u) + .9 (w) 17 (n ff ) ff2 M . (92) 
Applying the Langreth rules we find the lesser GF (omitting (w) for brevity): 

G< a = g< + g r U (n s ) g< + g<U (n s ) g a 2 (93) 
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where g < — g 7 "E,Qg a , g 2 = goods' an< ^ ^o" ^ s defined in equations and By definition g< is imaginary, 

hence, for G^ a to be imaginary one requires that 

Ax = g r U(n s )g^<g a 2 +g^<g a U(n s )g a 2 , (94) 

be imaginary. Since g 2 ^og 2 an d 9 r ^o9 a are P ure imaginary quantities, for G^ a to be pure imaginary, the real part 
of A\ needs to cancel, i.e.,: 

Im (g r ) g r 2 ^<g a 2 = -Im (g a 2 ) g r X<g a . (95) 

Define 

e = e a + Re(Y? Q ), (96) 
T = -Im (E5), (97) 

and use it to rewrite the equation (|95[) as: 

£ ^< \ _ zE. v]< - (98) 

(fiw - e ) 2 + (r) 2 ° (fiw - e - £/) 2 + (r) 2 (hiu - e a - U) 2 + (V) 2 ° (fcj - e ) 2 + (r) 2 ' 

Obviously the real part of A\ cancels, hence G^ CT is imaginary and fulfills the symmetry G< CT (uj) = — (G^ a (uj))* . 



VIII. FULL DERIVATION OF THE BROKEN SYMMETRY IN THE DOUBLE ANDERSON MODEL 



Now we refer to subsection The double Anderson model in the manuscript. In what follows we show in greater 
detail that (g^ (uj)) ^ ( {p^a. ) • Again, following the derivation in Ref. [5(3 we define the following contour 



ordered GFs: 



G% {t,t') = - l -(T c d atT (t)d\ a (t')), 



(99) 



-<t<x<x 

"a/37 



(t, t') = - % - (T c n aT (t) dp u (t) rft (*')) 



where t = a/a . The resulting EOM are: 



G% (uj) = (nw - - So Mr 1 x (5- + h a0 G% (uj) + U a G%% (uj) 
(w) + C/TG^ H) , 



(100) 



(101) 



[UJ 



[UJ 



[UJ 



[UJ 



(UJ) 



(tiu - e aa -U a - V£$ (n Pa ) - V$ (n fi9 ) - S (uj))' 1 

x [k p gi% (uj) + <rw) (v$g%% H + m)] , 

(ftj - ep a - ^ (n^) - V£ (n aa ) - V$ - S 

x [(n a9 ) + h° pa &Z a H + (n aff > (UpGppp (w) + V^G^ (w))] , 

(fej - £ p& - ^ W,) - V^f (n aff > - VJJ - S M)~ J 

x [<»«*) + fc£ G^ (w) + M (UpGfft (uj) + Vf°G%% (uj))] , 

(huj - e aa - U a (naff) - V$ (n Pa ) - V$ - E M) _1 

x 'hlpG 9 ^ (uj) + (npff) (u a G%% (uj) + VtfGffi (uj))] , 

(Jvjj - e a<7 - U a (n a ff) - V°p (rips) - V°p - E (uj)) 1 

x - (dl a d a , a ) + hlpGlffi (uj) + (n p>v ) (u a G%Z (uj) + V^Gf^ 

(huj - ep a -Up- Vf° (n as ) - V f % (n aa ) - S (uj))' 1 

M + h° Pa G%Z (uj) + (n ra ) fe ff G^ (uj) + V^G%^ (uj 



(uj) 



(102) 
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By applying the Langreth rules one can find the retarded and advanced projections of the single particle GF (equation 
(| 1 1 [) ) . For simplicity we derive them for the case where V? T = 0. Define (as usual omitting (ui) for brevity) 

{9i) r ' a = T / v r,a , (103) 

= fc 1 rj y r.a , (104) 

nw - Ei^ - Ui - L ' 



1 



and 



(105) 



where So is defined in equation (|50p. Now we are ready to look at the equation we get for ■ 

{G7pY = (9 a ) r h^(G^) r + {g a ) r U a {^) r , (106) 

(G^Y = (9pY + (9p) r hlAG^) r + (9pYU P (G$^) r , (107) 

(G^) r - (C) r ^(GS) r > (108) 

(G|^) r = (flfS) r ^/3(G|??)''+(flfS) r (n i9ff >Cr o (G^) r , (109) 

(G^) r = ( 5 0) r ( n/Jff ) + {g%) r h% a (G$^) r , (110) 

(K^Y =(5S?) r (^) r ^a(G^) r + (^) r (^)^(G|^) r . (Ill) 

Substituting (g^ into the equation of (g^\ : 

{G7pY = (g a ) r K 3 ((9pY + (gp) r (G^) r + ( 90 y u p {&$^) r ) + ( 9a ) r u a {gz^Y , 

(112) 

Now we need the equations for (^Z^j an< ^ {^pp^j • U sni g equations f| 108[) to (jllip we get: 

x (1 - (g$%) r h% a {g%l) T h^y 1 (g%l) r h% a {g%° a ) r {n [js ) U a ) ~* 

x ((gZY K,p {g%) r + (gZY K,p (gZpY M Up 

x (1 - {g$l) r h% a {g$l) r h^y 1 (g%) r (n 0s )) , (114) 



(Gf/y P Y = (!- {9%) T h°f>M«Y Kp- WY h U(9$*Y M U <* 

x (1 - (gZY K,p {g%Y h^y 1 ( 9 ZY K,p (aZpY (v) u p ) _1 

x {{g$lY M + {g$fY h Z (9%«Y M U a 

x (1 - {gZY K,p {9%Y h^y 1 (gZY K,p (gZ Y <na,*>) • (H5) 
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The same way one can derive an expression for I G aa 



with 



x ({wY h$ ia (g a ) a + (gp) a h% a (g a ) a U a (<G^) a + (gp) a Up (G%l) a ) , (116) 
(<S$7S = (l-(42)°^.a(ff|S) ^-(4?) ^,«(4S) <^>^a 

x (i - (gZ) a Kp {aUT h% a y X (aZT Kp (gZp) a M u p ) 1 

x {{MY h% a (g$l) a (n Pt9 ) + h% a {g%) a (np, 9 ) U a 

x {i-{gZ) a K {9Z) a hl a y\gZ) a {n aS )), (117) 



and 



a a.aa. ) 



= (i - (gZ) a Kp ( 9 Zp) a h% a - (gZ) a K {g%Y M Up 

x (i - ( g %) a hi a ( 9 $°) a hip)- 1 { 9 %y hi a ( g $iy ( n0 - a ) u a ) 1 

x ((gZT W + (gZT Kp {gZpY W u p 

x (i - {gfp) a hl a {g%yKp)~ x {9%%) a n, a {gfS (np,*)) . (118) 

The question we now ask is whether 

(GS5) r =((G£)°)\ (119) 

Using 

(«r =((*)")*, (flS) r =((fla) B )", (4) r =((4)T' ( i2 °) 

and looking back at equations (| 1 13|) and (| 1 we find that 

(/ - (g a ) r K,p (gp) r h^y 1 = ((/ - (gp) a hl a (g a ) a hl,p) - 1 ) * , (121) 

and that 

(g a ) r h^(gp j r = ((gp) a hl a (g a ) a y . (122) 
Therefore, it is sufficient check whether the next equality 

(g a ) r u a (GSS5) P + Cff«) r K,p (9pY Up {G$%) r = {{gp) a h^ a ( 9a ) a u a (g^)° + (gp) a Up > 

(123) 

holds. Substitute the equations for (g^Y , (Gf^Y , (G%£)" and fe^)" into equation (|T2"3) and after some 
tedious algebra we find that unless 

(4) r = (WT = W r , (124) 

the identity (G^j = (f^«) ) does not hold - But as 

^ = hu-ep, a -Up{np, s )-^ (125) 

fe) r = T p ~r, (126) 



and 



it is obvious that {g^p\ Y (i.9pY)* ■ hence finally \G%$J Y (O-'pa) ) ■ The same can be done to show that 
(G^) < ' > Y ((g?S) < " > ) and (G^) r - (G^)° ^ (G^) > - (G^) < 



